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ABSTRACT 

<f 12 

Observations of CO at high redshift indicate rapid metal enrichment in the nuclear regions of at least some galaxies 
. in the early universe. However, the enrichment may be limited to nuclei that are synthesized by short-lived massive 

y—i • stars, excluding classical "secondary" nuclei like 13 C. Testing this idea, we used the IRAM Interferometer to tentatively 

detect the 13 CO J=3— >2 line at a level of 0.3Jykms _1 toward the Cloverleaf Quasar at z = 2.5. This is the first 
observational evidence for 13 C at high redshift. The 12 CO/ 13 CO J=3— >-2 luminosity ratio is with 40lg 5 much higher 
than ratios observed in molecular clouds of the Milky Way and in the ultraluminous galaxy Arp220, but may be similar 
to that observed toward NGC 6240. Large Velocity Gradient models simulating seven 12 CO transitions and the 13 CO 
line yield 12 CO/ 13 CO abundance ratios in excess of 100. It is possible that the measured ratio is affected by a strong 
submillimeter radiation field, which reduces the contrast between the 13 CO line and the background. It is more likely, 
O . however, that the ratio is caused by a real deficiency of 13 CO. This is already apparent in local ultraluminous galaxies 

and may be even more severe in the Cloverleaf because of its young age (^2.5 Gyr). A potential conflict with optical 
data, indicating high abundances also for secondary nuclei in quasars of high redshift, may be settled if the bulk of the 
-£j CO emission is originating sufficiently far from the active galactic nucleus of the Cloverleaf. 
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et al. 2007; Juarez et al. 2009; Matsuoka et al. 2009). This 
evidence, mainly from optical lines, is supported by mil- 
limeter detections of CO and dust in high-redshift sources, At optical, near-IR, and UV wavelengths it is difficult 

indicating rapid metal enrichment due to starbursts in the to discriminate between an element's isotopes because their 

circumnuclear regions of at least some galaxies in the early atomic lines are blended (e.g., Levshakov et al. 2006). The 

universe (e.g., Solomon & Vanden Bout 2005). This en- prospects are better with radio lines from isotopic substi- 

richment, however, might apply mainly to atomic nuclei tutions in molecules, which are well separated by a few 

' that are synthesized in short-lived massive stars, and not percent of their rest frequency from the main species. This 

so much to "secondary" nuclei like 13 C that are thought separation allows both the main and rare species to be eas- 

to be mainly synthesized in longer-lived, less-massive stars ily identified, and to be observable with the same radio 

(but see, e.g., Hamann et al. 2002 for the mainly secondary receivers and spectrometers, 
element nitrogen). 

In the local universe, 12 C/ 13 C abundance ratios are _ , „ _ , , 

sometimes considered to be a diagnostic of deep stellar mix- The Cloverleaf Quasar (H1413+117), partly because of 

, c «_ • j) « j » i amplification by gravitational lensmg, is a mgn-z source 

mg and a measure of primary vs. secondary nuclear .,f J . ° . n , ... . To^e,^ n « 

processing (e.g., Wilson & Rood 1994). While 12 C is pro- with exceptional peak flux densities m C b O (hereafter 

duced by He burning on rapid time scales in massive stars, ^ see A PPendix 2 of Solomon & Vanden Bout 2005) 

13 n ■ -i t-u ■ a u • f i2n j 1 his source is therefore one of the best candidates to search 

3 C is mainly synthesized by CNO processing of C seed 13 16 13 

nuclei from earlier stellar generations. This processing oc- tor „, U , U l nereatter totT Y ™ te ft models ot cnem- 
ical evolution over a Hubble time. In this paper we report 

Send offprint requests to: C. Henkel, e-mail: on a search for 13 CO(3-2) emission in the Cloverleaf at 

chenkel@mpifr-bonn.mpg.de 2=2.5579, when the universe was 2.5 Gyr old. 
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Fig. 1. Contour map of continuum plus 13 CO J=3— >2 emis- 
sion, covering the central 960 kins -1 toward the Cloverleaf 
QSO. The beam is 6'/lx5'/4 (lower left) and the contour 
step is 0.09 mJy (ler). The peak value and the spatially- 
integrated intensity in the central source are both 0.6 mJy. 



2. Observations 

The measurements were made with the IRAM 
Interferometer on Plateau de Bure, France, in July, 
August, and September 2008, with 5 antennas in the 
compact D-configuration (maximum baseline 97 m) and 
the new dual-polarization receivers. The receiver and 
system single-sideband temperatures were 40 and 100 K, 
respectively. The spectrometers covered 1 GHz in each 
polarization, and the raw spectral resolution was 2.5 MHz, 
or 8.1 km s . The data were binned to various spectral 
resolutions; in this paper we present data binned in 
19x160 kms -1 channels, covering a range of 3040 kms -1 , 
with a noise of 0.22 mJy beam -1 (ler) in each channel. The 
naturally-weighted synthesized beam was 5 / .'6x4'. / 8 at p. a. 
62°. Because the four CO spots of the lensed Cloverleaf 
image are spread over l'/7, we included more of the total 
flux by applying to the u, v data a Gaussian taper that fell 
to 1/e at a radius of 100 m. The slightly broadened beam 
then became 6'.'lx5'.'4, and the noise in the individual 
channels is 0.23 mJy beam -1 . 

3. Results 

Figures 1 through 3 show the data, and Table 1 summa- 
rizes the results. In the integrated line -I- continuum map 
(Fig. 1), the peak position (Table 1) agrees well with the 
centroid of previous high-resolution interferometer maps of 
the source (e.g., Alloin et al. 1997; Yun et al. 1997; Kneib 
et al. 1998). At 93 GHz, the expected continuum is 0.30- 
0.35 mJy (from Fig. 3 of Weifi et al. 2003 and the power-law 
given in Bradford et al. 2009) and a map in the 13 off- 
line channels at the positive and negative velocity ends of 
our spectra indeed yields a continuum flux of 0.3±0.1mJy 
(Fig. 2). This continuum adds to the line signal, and for 
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Fig. 2. Contour map of the 3.2 mm continuum emission, 
covering 2080 kms -1 in the off-line channels. The beam 
is 6'.'lx5'.'4 (lower left) and the contour step is 0.09 mJy 
(1.5 a). The peak and the spatially-integrated flux den- 
sity of the central source are both 0.3 mJy. Combining 
this figure with the previous one, we conclude that in the 
960 kms -1 band centered on the 13 CO J=3— >2 line, line 
and continuum each contribute half of the total flux density. 



this reason, the line appears much broader than the ~ 430 
kms -1 widths of the 12 CO and [Ci] lines (Weifi et al. 2003). 
The observed line has a low signal-to-noise ratio, which pre- 
vents a clear distinction between line and continuum, and 
does not allow us to constrain the line shape. Above the 
0.3 mJy continuum, a Gaussian fit yields an integrated line 
flux of (0.3±0.1) Jykms -1 (Fig. 3, see also the much higher 
upper limit given by Barvainis et al. 1997, their Table 1). 
An alternative Gaussian fit, with the line width fixed to the 
width of the 12 CO line, yields a peak line flux density of 
(0.44±0.12) mJy beam -1 , and the same integrated line flux 
as the fit shown in Fig. 3. This integrated flux, corrected 
for frequency squared, leads to a 12 CO/ 13 CO J=3— >2 line 
luminosity ratio (= brightness temperature ratio) of 40^g 5 
(Table 1). This value is conservative. With the line width 
fixed to the width of the 12 CO line and the actual peak flux 
density of order 0.35 mJy, the ratio would become ^75. 

4. Large velocity gradient model calculations 

12 CO lines have higher optical depths than those of 13 CO. 
Therefore, the measured 12 CO/ 13 CO line intensity ratio 
(Sect. 3) is a lower limit to the 12 CO/ 13 CO abundance ra- 
tio. To further constrain the 12 CO/ 13 CO abundance ratio 
of the Cloverleaf QSO, Table 2 provides flux densities and 
brightness temperatures of seven 12 CO transitions. To sim- 
ulate these values, a Large Velocity gradient (LVG) model 
was used with collision rates from Flower (2001), a cos- 
mic microwave background of 9.7 K, and an ortho-to-para 
H2 abundance ratio of three (e.g., Weifi et al. 2005, 2007; 
Riechers et al. 2006b). The latter is, however, not critical 
for this study. 
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Fig. 3. CO J=3^2 from the Cloverleaf QSO, measured with the IRAM interferometer. Left: 12 C0(3-2) profile in 
10 kms' 1 channels from Weiss et al. (2003). Velocity offsets are relative to 97.1928GHz Right: 13 CO(3-2) profile in 
160 kms -1 channels, from this paper. Velocity offsets are relative to 92. 91816GHz (z — 2.55784). The red curves show 
Gaussian fits above a continuum of 0.3 mJy. 



Table 1. 13 CO(3-2) Observations and results. 



Table 2. CO line ratios in the Cloverleaf. 



Parameter 



Observed CO(3-2) quantities: 

R.A. (J2000) 

Dec. (J2000) 

Center frequency (GHz) 

Redshift (LSR) a) 

Continuum flux density (mJy) 

Integrated 13 CO flux (Jy kms" 1 ) 

Derived CO (3-2) quantities: 
L'( 13 CO (K kms-'pc 2 ) 1 ' 
Z/( 12 CO (K kms-'pc 2 ) 1 ' 
L' ratio 12 CQ/ 13 C0(3-2) 



j CQ(3-2) 



Line 



14 h 15 m 46=28 ± 0=03 
+ 11° 29' 44'.'0 ± 0"4 

92.91816 
2.55784 ± 0.00003 

0.3 ±0.1 

0.3 ±0.1 



(1.1 ±0.3)xl0 10 
(45.9±3)xl0 10 



40 



+ 25 



a) adopted from 12 CO (Weiss et al. 2003). 

b) in a beam of 6'.'lx5'.'4. 

°' This is the lens-amplified value for a luminosity distance 
of D L = 21.28 Gpc (H =71 kms^Mpc -1 , Q, m = 0.27, 
fivac = 0.73) and an angular diameter distance of Da = 
1.682 Gpc; linear scale: 1" O 8152 pc (Wright 2006). 



CO (3-2) 
CO (4-3) 
CO (5-4) 
CO (6-5) 
CO (7-6) 
CO (8-7) 
CO (9-8) 



Integrated 
line flux 
(Jy kms -1 ) 



T b ratio a 
to 

12 CO(3-2) 



Reference ' 



13.2±2.0 
21.1±3.2 
24.0±3.6 
37.0±5.6 
45.3±6.8 
51.4±7.7 
41.8±6.3 



1.00±0.15 
0.90±0.13 
0.65±0.09 
0.70±0.10 
0.63±0.09 
0.55±0.08 
0.35±0.05 



3 CO(3-2) 0.3±0.1 



0.025 



-0.000 



a ' If all lines have the same area filling factor. Adopted 
la errors are ±10% for the flux densities and ±15% 
for the brightness temperature ratios. 

b) (1) Weifi et al. (2003); (2) Barvainis et al. (1997); 
(3) Bradford et al. (2009); (4) this paper. 



We calculated a grid for 12 CO/ 13 CO with kinetic 
temperatures of 30-100 K and 12 CO fractional abun- 
dances per velocity interval of [ 12 CO]/([H 2 ](dv/dr)) = 
10~ 4 '~ 7 pc (kms -1 ) -1 . Accounting for possible effects of 
cloud structure, not only a spherical but also a plan- 
parallel cloud morphology was considered, with escape 
probabilities /3 sph cricai = (l-e~ r )/r and /3 p i a n-paraiici = 
(l-3e -r )/(3T), respectively (t: optical depth). Resulting 
12 CO/ 13 CO abundance ratios reproducing the six measured 
12 CO line intensity ratios (Table 2) are given in Figs. 4 and 
5 together with reduced x 2 (x 2 e d) va l ues of the best fit. We 
adopted a la error of 15% for each fitted brightness tem- 



perature ratio. The dependence of the resulting 12 CO/ 13 CO 
ratios on cloud morphology is caused by the different escape 
probabilities, related to r in the case of a spherical and to 
3r in the case of a plan-parallel cloud geometry. Therefore, 
a required amount of excitation through photon trapping 
is reached at lower 12 CO optical depths in the case of a 
plan-parallel morphology, resulting in smaller 12 CO/ 13 CO 
abundance ratios. 

The Xred va l ues displayed in Figs. 4 and 5 indicate that 
the CO data can be fitted by a single molecular gas com- 
ponent (cf. Bradford et al. 2009). All calculations are also 
consistent with the (not very stringent) upper limits for the 
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13 CO J = 7— >6 and 8-»7 flux densities from Bradford et 
al. (2009). At first sight, the figures do not strongly reduce 
the permitted parameter space, providing x 2 ed values of 
order 1.25-2. In the upper left corners of each figure, how- 
ever, the Xred va l ues rise significantly, becoming too large 
to provide credible solutions. As a consequence, the over- 
all 12 CO/ 13 CO abundance ratio appears to be >100 in the 
Cloverleaf QSO. There exist further constraints: (1) Tkin < 
30 K is prohibitive because of the temperatures determined 
from Ci (-30 K) and the dust (-50 K, Weifi et al. 2003). 
Furthermore, such low temperatures would require extreme 
CO column densities to raise photon trapping to such lev- 
els that the emission from the higher J transitions could 
be reproduced. (2) Tkin > 50 K is also not likely because of 
the temperature deduced from [Ci] and the close associa- 
tion of CO and Ci, which appears to be independent of the 
environment (e.g., Ikeda et al. 2002; Zhang et al. 2007). 

For 30<T kin <50K and [ 12 CO]/([H 2 ](dv/dr)) = 
10~ 7 pc (kms -1 ) -1 , we obtain 12 CO/ 13 CO abundance 
ratios in the range 200-3000 (Figs. 4 and 5). However, 
such a low fractional abundance per velocity interval 
can be firmly excluded. With [Ci]/[H2] reaching values 
in agreement with those of the local galactic disk (Weifi 
et al. 2005), the [ 12 CO]/[H 2 ] abundance ratio should be 
of order 10~ 4 (e.g., Frerking et al. 1982). The resulting 
velocity gradient of dv/dr = 10 3 kms^ 1 pc -1 would be far 
too large in view of the measured line width (e.g., Weifi 
et al. 2003) and the kinetic temperature such extreme 
conditions would induce (e.g., Wiklind & Henkel 2001). 
A velocity gradient of a few km _1 pc _1 is more realistic 
as, e.g., obtained from clouds in virial equilibrium for 
densities of order 10 4 cm~ 3 (Goldsmith 2001; his equation 
2). Such densities are commonly derived for high z sources 
(e.g., Weifi et al. 2005, 2007). For diffuse clouds, velocity 
gradients should be larger (e.g., Papadopoulos et al. 2010). 
Bradford et al. (2009) suggest that in the Cloverleaf the 
velocity dispersion may exceed the virial requirement by at 
least an order of magnitude. Therefore the best choice may 
be [ 12 CO]/([H 2 ](dw/dr)) = lO" 5 - 6 pc (kms" 1 )" 1 (for 
the higher value see, e.g., Riechers et al. 2006b; Weifi et 
al. 2007) to simultaneously fit the observed CO transitions 
from J=1^0 up to 11—^10. Depending on the adopted 
kinetic temperature (30-50 K) and cloud morphology, and 
irrespective of the optimal [ 12 CO]/([H 2 ](du/dr)) value (as 
long as it is in the wide range displayed by Figs. [4] and [5]) 
we then find a 12 CO/ 13 CO abundance ratio in the range 
300-10000. In the following we will discuss whether this 
estimate can be realistic. 



5. Discussion 

In order to further evaluate our observational result, 
we have to discuss the correlation between molecular 
12 CO/ 13 CO and atomic 12 C/ 13 C abundance ratios and to 
summarize relevant observational data from low-redshift 
galaxies which are, like the Cloverleaf, ultraluminous in the 
infrared. Finally, we will address some fundamental prob- 
lems, which are related to the still poorly known morphol- 
ogy of the gas surrounding the Cloverleaf QSO. 
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Fig. 4. Results from large velocity gradient (LVG) radia- 
tive transfer calculations using a spherical cloud model to 
simulate the line intensity ratios given in Table 2. The com- 
mon logarithm of the 12 CO/ 13 CO abundance ratio is shown 
as a function of kinetic temperature in units of Kelvin 
and of fractional abundance in units of pc (kms -1 ) -1 for a 
12 CO/ 13 CO J=3->2 line intensity ratio of 40. Resulting re- 
duced x 2 values (x 2 d) f° r tne simulation of the four 12 CO 
lines given in Table 2 are shaded. Lightest grey: 1.00 < x? e d 
< 1.25, darker shades of grey at 1.25, 1.50, 1.75 to 4.00 with 
an increment of 0.25. The maximum value in the upper left 
corner is x 2 e( j = 4.07. 

5.1. Chemical fractionation and isotope selective 
photodissociation 

Observed isotope ratios may be affected by fractionation. 
The 12 CO/ 13 CO abundance ratio is likely influenced by the 
reaction 

i3 C++ i2 C0 ^ 12 C+ + 13 CO + A^ 35K 

(Watson et al. 1976). The process enhances 13 CO relative to 
12 CO in the more diffuse C + rich parts of molecular clouds. 
This may be compensated by isotope selective photodis- 
sociation. 12 CO and 13 CO need similar amounts of self- 
shielding to survive in a hostile interstellar environment. 
This favors the more abundant isotopologue (e.g., Sheffer 
et al. 2007). 

For the Galaxy, such effects can be quantified. Milam 
et al. (2005) summarized 12 C/ 13 C ratios from the galac- 
tic disk, obtained with the three molecules CO, CN, and 
H 2 CO. These molecular species are synthesized by quite 
different chemical reactions. The good agreement between 
their 12 C/ 13 C ratios and a lack of correlation with kinetic 
temperature suggests that chemical fractionation as well as 
isotope selective photodissociation do not greatly affect the 
determined isotope ratios. 

Whether this result is also valid in the case of the 
Cloverleaf QSO may not be obvious at first sight. The ul- 
traviolet radiation field in the vicinity of the quasar might 



Henkel et al.: CO in the Cloverleaf Quasar 



5 




7 log [([ 12 CO]/[H 2 ]) * (dv/dr)" 1 ] * 



Fig. 5. Same as Fig. 4, but for a plan-parallel cloud ge- 
ometry. Resulting x 2 ec i va l ues f° r the simulation of the line 
intensity ratios given in Table 2 are shaded. Lightest grey: 
1.00 < x 2 < 1.25, darker shades of grey at 1.25, 1.50, 1.75 
to 4.00 with an increment of 0.25. The maximum value in 
the upper left corner is x 2 ec i = 6.31. 

be exceptionally strong, favoring 12 CO over 13 CO and thus 
leading to an enhanced molecular abundance ratio with re- 
spect to 12 C/ 13 C. However, such a scenario is not likely. 
Firstly, most of the galactic data were obtained toward 
prominent sites of massive star formation, where the UV 
radiation field is also exceptionally intense. Secondly, judg- 
ing from Cl, in the Cloverleaf the excitation of the molecular 
gas is intermediate between conditions found for the star- 
burst galaxy M82 (T ex ci ~50K) and the central region of 
the Milky Way (T cx , C i ~ 22 K) (Stutzki et al. 1997; Weifi et 
al. 2003). Thirdly, polycyclic aromatic hydrocarbon (PAH) 
features are as strong as expected with respect to the far in- 
frared luminosity when compared with more nearby ultralu- 
minous star-forming galaxies, favoring "normal" conditions 
and a predominantly starburst nature of the Cloverleaf's 
huge FIR emission (Lutz et al. 2007). Finally, the CO emis- 
sion from the Cloverleaf appears to be more extended than 
the effective radius out to which the quasar could dominate 
the UV field. 

Modeling both the source and the lens of the Cloverleaf 
QSO, Venturini & Solomon (2003) find a characteristic 
radius of r ~ 800 pc for the CO J=7-6 line, which is 
higher excited and thus possibly less widespread than the 
J=3— >2 transition considered here. If the Cloverleaf's in- 
trinsic far infrared luminosity (Lfir ~ 5xl0 12 L Q , Lutz 
et al. 2007) would entirely originate from 6.2-13.6 eV pho- 
tons emitted by the active nucleus, we would obtain, at 
a radius of 800 pc, a UV photon illumination of x ~ 
10 5 xo with respect to the local galactic radiation field, 
Xo = 2xl0~ 4 crgcm" 2 s" 1 sr" 1 (see Draine 1978). The 
Cloverleaf QSO is a Broad Absorption Line (BAL) quasar 
which permits at least a partial view onto its nuclear en- 
gine. Therefore, taking the Cloverleaf's UV luminosity from 



Fig. 1 of Barvainis et al. (1995) and accounting for a gravi- 
tational amplification by a factor of 11 (Solomon & Vanden 
Bout 2005), we obtain accordingly x ~ 2.5xl0 4 xo- Both x 
values are consistent with those encountered in prominent 
galactic sites of massive star formation and may be up- 
per limits if the Cloverleaf posseses a self-shielding rotating 
disk. To summarize, physical conditions in the Cloverleaf 
host galaxy appear to be sufficiently normal so that the 
12 C/ 13 C isotope ratio should not strongly deviate from the 
12 CO/ 13 CO molecular abundance ratio. 

5.2. 12 CO/ 13 CO ratios in z<l galaxies 

In our Galaxy, the 12 CO/ 13 CO line intensity ratios from 
molecular clouds are typically about 5, probably corre- 
sponding to true 12 C/ 13 C abundance ratios of ^25 in the 
galactic Center, ~50 in the inner galactic disk and the 
LMC, ^70 at the Sun's galactocentric radius, and ^100 
in the outer Galaxy. The solar system ratio of 89 may have 
been typical of the galactic disk at the Sun's galactocentric 
radius 4.6 Gyr ago (e.g., Wilson & Rood 1994; Wouterloot 
& Brand 1996; Wang et al. 2009). Within the framework 
of "biased infall", where the galactic disk developed from 
inside out (Chiappini et al. 2001), there might be a future 
chance to use 12 C/ 13 C ratios as a chronometer for nucle- 
osynthesis. 

In nearby galaxies, the 12 CO/ 13 CO line intensity ratios 
are usually measured in the J=l-0 line and have typical 
values of ^10. They are higher than the values for indi- 
vidual molecular clouds in the Galaxy because they are 
mostly observed with larger beams. These include not only 
the dense clouds, where both species are (almost) optically 
thick, but also the molecular intercloud medium, where 
13 CO is optically thin. Like the better-resolved CO line 
ratios in our Galaxy, the ratios in nearby galaxies probably 
correspond to true 12 C/ 13 C abundance ratios between 40 
and 90 (e.g., Henkel et al. 1993). 

In a presumably "normal" spiral galaxy at redshift 0.89, 
in the lens of the background source PKS 1830-211, Wiklind 
& Combes (1998), Menten et al. (1999), and Muller et al. 
(2006) derive, from the optically thin wings of the absorp- 
tion lines of HCO+, HCN, and HNC, a 12 C/ 13 C abundance 
ratio of 27±2. Apparently, even at an age of the universe of 
~6.5Gyr, it appears that 13 C is as abundant with respect 
to 12 C as in the center of our Galaxy at the present epoch. 

Some low-redshift (ultra)luminous infrared galaxies 
((U)LIRGs), however, show peculiarities, which may be rel- 
evant to the Cloverleaf. Local (U)LIRGs are known to re- 
veal 12 CO/ 13 CO J= 1— >0 line intensity ratios which tend 
to be higher than the canonical value of 10 for "normal" 
galaxies (see, e.g., Aalto et al. 1991; Casoli et al. 1992; 
Henkel & Mauersberger 1993). According to Taniguchi & 
Ohyama (1998), there is a tight correlation between L( 12 CO 
J=l— >0) and £fir- However, when comparing "normal" 
galaxies with those with a high 12 CO/ 13 CO J=l— >0 ra- 
tio, the 13 CO luminosities show a deficiency by an average 
factor of ^3, This 13 CO deficiency is readily explained by 
metallicity gradients in the progenitor galaxies and strong 
interaction- or merger- induced inflow of gas into the lu- 
minous cores (e.g., Rupke et al. 2008). Apparently, for ul- 
traluminous galaxies the common luminosity - metallicity 
correlation is not valid. Ultraluminous galaxies are charac- 
terized by a lower metallicity, likely yielding higher 12 C/ C 
isotope ratios. In the early universe, gas from outside the 
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cores of the merging progenitors may have been particularly 
metal poor, leading to extreme carbon isotope ratios. 

For T kin ^ 20 K, the 12 CO J=3-^2 line is more opaque, 
typically by a factor of 3, than the corresponding l—>0 line. 
Thus our conservatively estimated J=3-2 12 CO/ 13 CO line 
intensity ratio of ^40 jig 5 corresponds to a 1—5-0 ratio well in 
excess of 40. So far, only few 12 CO/ 13 CO J=3->-2 line ratios 
have been measured in luminous mergers of low redshift. 
Greve et al. (2009) find 8±2 for the ULIRG Arp 220 and 
>30 for the LIRG NGC 6240. The latter value might be 
consistent with that of the Cloverleaf. 

5.3. Are there alternatives to a 13 C deficiency in the 
Cloverleaf? 

Sects. 5.1 and 5.2 suggest, that our measured 12 CO/ 13 CO 
line intensity ratio (or its lower limit) require a significant 
13 C deficiency in the Cloverleaf. Are there caveats we may 
have overlooked when reaching this conclusion? 

If the bulk of the CO emission would not arise, as 
suggested by Venturini & Solomon (2003), from a molec- 
ular disk but from a large scale outflow, such gas would 
not be in virial equilibrium and could arise predominantly 
from a diffuse gas phase. While this would yield (within 
the LVG approach) a higher velocity gradient and a lower 
[ 12 CO]/([H2](dv/dr)) value than what is needed for viri- 
alized clouds, required densities would then be well in ex- 
cess of 10 4 cm~ 3 , in contradiction with our assumption of 
predominantly diffuse gas. Furthermore, as long as Tki n re- 
mains moderate (^50 K; see Figs. S] and [5]), 12 C/ 13 C ratios 
remain larger than those encountered in the galactic disk 
(Sect. 5.1). 

Following White (1977), radiative transfer models with 
simple geometry, either based on microturbulence or on 
systematic motions, lead to peak and integrated intensi- 
ties which agree within the differences (up to a factor of 
three) caused by an uncertain cloud geometry. A full 3-D 
model of a rotating circumnuclear disk, computing the ra- 
diative transfer through many lines of sight, calculating the 
LVG level populations within each pixel of the simulated 
source, and also including continuum radiation from dust 
(e.g., Downes & Solomon 1998) may be worth doing. In the 
Cloverleaf, however, the distribution of the molecular gas 
is still poorly known. 

A large 12 C/ 13 C ratio, implying an underabundance of 
13 C, appears to be in direct conflict with optical data. As 
already mentioned in Sect. 1, solar or super-solar metallici- 
ties are common in quasars up to high redshifts. This does 
not only refer to so-called "a-elements" being rapidly syn- 
thesized in short-lived massive stars but also to iron (e.g., 
Iwamuro et al. 2004; Kurk et al. 2007; Sameshima et al. 
2009), carbon (e.g., Jiang et al. 2007; Juarez et al. 2009), 
and, even more importantly, nitrogen (Hamann & Ferland 
1999; De Breuck et al. 2000; Vernet et al. 2001; Hamann 
et al. 2002; Nagao et al. 2006; Matsuoka et al. 2009), with 
14 N being mainly a secondary nucleus produced by CNO 
burning just like 13 C. A possible explanation for the con- 
tradictory results obtained at or near optical wavelengths 
and the microwave data presented here may be different lo- 
cations. It is well possible that mainly secondary nuclei like 
13 C and N are enriched close to the quasar, in the Broad 
and Narrow Line Regions and in outflows originating from 
the active galactic nucleus (AGN). However, CO J=7— >6 



may arise hundreds of pc away from the AGN (Venturini 
& Solomon 2003) and some of the J=3— >2 photons may be 
emitted from locations even farther away. 

There exists, however, also the possibility that our mea- 
sured high 12 CO/ 13 CO luminosity ratio is misleading and 
does not imply a large 12 C/ 13 C ratio. As a consequence 
of different optical depths, 12 CO lines are almost thermal- 
ized and are characterized by excitation temperatures well 
above the level of the cosmic microwave background even 
at z=2.5. 13 CO is less thermalized. In our best fitting mod- 
els, its J=3— >2 excitation temperature lies in the range 
20-30 K. This is significantly above the 9.7 K of the CMB. 
However, an extreme (and therefore unlikely) enhancement 
of the background level by dust radiation could reduce the 
contrast between line and background for 13 CO far more 
efficiently than for 12 CO (see Papadopoulos et al. 2010 for 
the case of Arp 220), thus establishing an apparent 13 CO 
deficiency. 

6. Outlook 

Molecular lines from galaxies in the distant universe have 
the potential to reveal the contribution of early stellar gen- 
erations to the enrichment of the interstellar medium. Our 
data from the z — 2.5 Cloverleaf QSO are a first step 
toward studying the isotopic composition of such gas in 
the distant past. Our data indicate, not unexpectedly, a 
strong deficiency of 13 C with respect to 12 C in the host 
galaxy. However, the weakness of the tentatively detected 
line, the limited number of observed transitions, the poorly 
constrained source morphology, and the potential influence 
of an enhanced submillimctcr radiation background do not 
yet allow us to derive a definite 12 C/ 13 C isotope ratio. 
Significant progress in this field either requires the detection 
of stronger sources or the higher instrumental sensitivity of 
the Atacama Large Millimeter Array (ALMA), which will 
allow us to study the isotopes of C, N, and O in a number of 
highly rcdshifted targets. Toward the Cloverleaf, the main 
isotopologes of HCN, HCO+, and CN (Solomon et al. 2003; 
Ricchers et al. 2006a, 2007) have already been detected. 
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